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PREFACE 

In  the  present  supplement  to  the  Manual  on  Supersonic 
Plow  and  Shock  Waves   the  gas  dynamical  phenomena  of  one- 
dlmenslonal  flow  Involving  detonation  and  combustion  are 
analyzed  from  the  mathematician's  viewpoint.     As   In  the 
Manual, content  and  emphasis   In  the   present  supplement  are 
conditioned  by  the  background  from  which  the  writers   happened 
to  approach  the   subject J      Important  points,   such  as   the 
finite  width  of   the  reaction  zone,   are  touched  upon  only  In 
an  appendix  and  In  the  bibliography. 

Jouguet,   In  his   classical  work,  was   concerned  mainly 
with  the  discussion  of  the  discontinuous  reaction  front; 
but  only  the   consideration  of  the  flow  as  a  whole  can  supply 
the   Information  necessary  to  determine   the  dynamical  phenomena 
Involving  detonation  or   combustion.     G.  I.    Taylor  has   studied 
detonation  processes  under  such  aspects.     In   the  recent 
research  program  of   the  Applied  Mathematics  Group  of  New  York 
University  a   somewhat  more   systematic  analysis   of  the 
mathematical  possibilities  of  flows   Involving  reaction  procesaei 
became  desirable,   and  the  present  report  gives  an  account  of 
such  investigations. 

R.   Courant 
Director,   Contract  OEMsr-945 
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REMARKS  ON  THE  MATHEMATICAL  THEORY  OF  DETONATION 
AND  DEFLAGRATION  WAVES  IN  GASES 

(Supplement  to  the  Manual  on  Supersonic 
Flow  and  Shock  Waves ) 

AMP  Report  38. 3R 


In  a  process  of  detonation  as  well  as  of  combustion 
(or  deflagration)  a  "reaction  front"  sweeps  over  the  com- 
bustible substance,  here  assumed  as  a  gas,  and  separates 
the  explosive  from  the  bxirnt  gases.   As  seen  In  the  "Manual 
on  Supersonic  Flow  and  Shock  Waves"  (AMP  Report  38. 2R)  there 
Is  a  close  analogy  between  "shock  fronts"  and  "reaction 
fronts,"  both  are  discontinuity  stirfaces  across  which  the 
state  of  the  substance  undergoes  sudden  changes  by  the  three 
laws  of  mass,  of  momentum,  and  of  energy;  the  difference  Is 
merely  that  across  reaction  fronts  the  energy -balance  con- 
tains as  one  term  the  energy  liberated  by  the  reaction,  while 
no  such  term  appears  In  the  energy  condition  for  shocks.* 

Together  with  the  differential  equations  In  the  zones 
of  continuity  and  with  Initial  and  boundary  conditions,  the 
transition  conditions  should  determine  the  flow  altogether. 
However,  while  such  a  determinancy  exists  for  gas  motions 
Involving  shocks  only,  a  different  situation  presents  itself 
as  soon  as  "reaction  fronts"  occur  (we  shall  use  this  term 
Indiscriminately  for  the  transition  front  due  to  detonation 
or  deflagration);   then  the  conditions  as  usually  envisaged 
In  non-reactive  gas  flows  are  no  longer  svifflclent  to  determine 
the  flow,  and  additional  conditions  must  be  found;  In  the  case 
of  deflagration  or  combustion  the  degree  of  indeterminacy  is 
even  higher  than  for  detonation  fronts. 

»  See  Manual  [l].  Chapter  III,  Section  50. 
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Thus  the  following  general  problem  arises:  To  what 
extent  are  reaction  waves  left  undetermined  by  the  laws  of 
conservation?  ?rhat  are  appropriate  additional  conditions 
to  determine  the  reaction  front?  These  questions  are  of  a 
strictly  mathematical  character.   They  refer  to  the  determinacy 
of  the  solutions  of  the  underlying  differential  equations  satis- 
fying well-defined  transition,  initial,  and  boTindary  conditions. 

Iii5)licitly  the  assumption  seems  to  have  always  been  made 
that  for  detonation  waves  the  well  known  Chapman- Jouguet  rule 
as  an  additional  condition  suffices  to  determine  the  process 
while  for  combustion  fronts  other  conditions  have  to  be  supplied 
from  a  specific  theory  of  the  mechanism  of  combustion.  Our 
following  mathematical  analysis  has  the  objective  of  clarifying 
this  statement  and  of  specifying  which  reaction  fronts  are  com- 
patible with  conditions  as  they  actually  may  occur.   Moreover, 
it  will  be  seen  in  what  manner  reaction  fronts  can  and  must  be 
accompanied  by  rarefaction  waves,  compression  waves  or  shocks. 
The  mathematical  analysis,  restricted  here  to  the  case  of  one- 
dimensional  motion  in  gases,  will  be  based  on  the  theory  of 
characteristics;  in  particular,  on  the  concept  of  domain  of 
dependence* 

While  the  following  discussion  stresses  the  mathematical 
aspect  of  the  theory,  it  may  have  physical  significance  in  parti- 
cular for  the  understanding  of  reaction  processes  taking  place 
in  gases  which  are  in  non -uniform  motion. 

Remarks  on  the  mathematical  theory  of  detonation  and  deflagraticn 
waves  in  gases* 

1.  Basic  notions  concerning  reaction  fronts. 

We  first  recall  a  number  of  well-known  facts  concerning 
reaction  fronts  in  gases.   The  gases,  before  and  after  reaction 
are  assumed  ideal  but  not  necessarily  polytropic.  In  other  words, 

♦   See  Courant-Hllbert  [2],  Chapter  5,  and  the  Manual  [1],  Chapt II. 

**  For  the  physico-chemical  notions  and  facts  in  the  present  re- 
port Lewis  and  v.  Elbe  [3].  For  the  theory  of  reaction  fronts 
see  Jouguet  [4] ,[5],  Becker  [6],  v.  Neumann  [7], Taylor  [14]. 

CONFIDENTIAL 


CONFIDENTIAL 


the  internal  energy  Is  a  function  of  the  temperature  but  not 
necessarily  proportional  to  lt»   By  p,  />,  r  =  p      ,    c  we 
denote  pressure,  density,  specific  volume,  and  sound  speed  of 
the  burnt  gas  quantities  referring  to  the  unburnt  gas  or 
"explosive"  are  characterized  by  a  subscript  (o).   We  Intro- 
duce the  quantity 

9  =rp, 

which  Is  proportional  to  the  temperature  T.  By  e,  1  =  e  +  0, 
g  ne  denote  respectively  the  Internal  energy,  enthalpy,  and 
energy  of  formation   per  unit  mass.   The  quantities 

E  =  e  +  g,   I=l  +  g 

may  conveniently  be  called  "chemical  energy"  and  "chemical 
enthalpy"  p.u.  mass  respectively.  We  ass\ime  that  these  quanti- 
ties depend  only  on  the  temperature  and  hence  "  on  9. 

2 
The  sound  speed  c  is  defined  by  the  relation  dp  -c  dyo= 

0  if  dE  +  pdr  =  0;  the  latter  relation  expresses  that  p  and  t 
are  coupled  adlabatlcally.  We  introduce  a  function  i=    ^(9) 
of  9  by  the  relation 


(1) 


dE  _   1 
35  -  7^ 


»   9  =  m'T^T,  R  being  the  gas  constant,  M  the  molecular  weight. 
Since  the  composition  of  the  burnt  gas  and  hence  its  molecu- 
lar weight  depends  on  the  temperature,  the  dependence  of  9 
on  T  is  not  strictly  linear. 

**  With  reference  to  absolute  zero  temperature.   By  liberated 
energy  f  we  shall  denote  the  difference  go-g  of  the  energies 
of  formation  at  absolute  zero  temperature  (not,  as  more 
customary,  at  a  standard  temperature  such  as  298°K. ) 

***  Strictly  speaking  this  is  not  correct  for  the  burnt  gas, which 
at  different  ten^jeratures  has  a  different  chemical  composi- 
tion, since  the  state  of  equilibrium  depends  not  only  on  the 
temperature  but  also  to  a  certain  degree  on  the  pressure.   In 
the  following  we  ignore  this  dependence. 
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Then  (see  App.  II)  the  sound  speed  is  given  by^ 

(2)  c^  =^©   . 

We  assume  the   gas   flow   to  be   one -dimensional,   taking 
place   in  a   long  cylindrical   tube.      By  u,  ^ ,    and 

(3)  V  =  u  -   ^ 

we  denote  gas  velocity,  reaction  front  velocity,  and  relative 
velocity  respectively. 

The  laws  of  conservation  in  the  regions  where  the  flow 
is  continuous  and  adiabatic  are  expressed  by  the  differential 
equations 


E^.  +  uE^  +  p{T^  +  ur^) 


Across  a  discontinuity  front  the  laws  of  conservation  are 

I  P-"     =/'o%  * 

II  P  ^p^^   =  Po  ^  Po^%      > 
m                                      |v2  -.I=|v2^I^  . 


*  For  we  have  identically 

dE  +  pdr  =  -i—  (rdp  +  )?pdr)  =  :T~-(dp  -  ^©d^), 
3r  -  1  J  -  1 

hence  vanishing  of  dE  +  pdr  Inplies   dp/d/o    =^Q. 
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The  mechanical  relations  I  and  II  Imply  relations 


(4)  '^^-A'-fK 

'o 


p-p 


(5) =  vv 


(6)  (p-Po)(  ^o-^^  =  <^-^o^^  ' 


(7)  o-  =  ^v  =  /o^v^ 


to  which  we  shall  refer  frequently. 

To  analyze  the  significance  of  the  thermodynamlcal  condi- 
tion III  we  eliminate  v  and  v^  and  obtain 

IV'  I^®^  -  Iq  =  i  ^  '^  *  '^o^  ^P  ■  Po^  °^ 


IV  E^®^  -  ^0=  "I  ^'^  '   '^o^  ^P  ■*■  Po^ 


Equation  IV  represents  a  relation  between  r    and  p,  when 
/5q,  p^,  and  hence  9^  are  kept  fixed.   Its  graph  in  a  (  r,p) 
-plane  is  indicated  in  the  diagram.   For  polytropic  gases  with 
y  =  const.,  the  graph  is  evidently  a  hyperbola,  since  then  E 
is  proportional  to  pT  ,  as  in  I. 

As  will  be  seen  later,  (see  Fig. 4),  it  is  only  to  the 
branches  r<  r^   and  to  the  branch  p  -  Pq  that  real  values  of  the 
velocities  v  ,  v  correspond.   Every  point  on  these  two  branches 
corresponds  to  a  possible  state  of  the  burnt  gas  after  it  has 
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*} 
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>s^ 

Ko 
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1 
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Figure  1 

Relation  between  pressure  and  specific 

volume  behind  a  reaction  front. 


undergone  a  transition  through  a  reaction  front  from  a  state 
(Tq,P^)  which  characterizes  the  unexploded  gas.  A  transition 
to  a  state  characterized  by  a  point  on  the  branch  t5  t^  is 
customarily  identified  with  a  detonation.   The  branch  p  -  p^ 
corresponds  to  deflagration  or  combustion.  The  transition  to 
P^,  '^ p^   =  T  ,  E^  =  E  corresponds  to  a  "constant  volume  detona- 
tion"; (in  this  limiting  case  the  velocity  v^  is  infinite)  the 
transition  to  Pr,  =  p  ,  f   i  =  i   corresponds  to  a  constant 


B^ 


pressure  combustion  (here  the  velocities  v  and  v  are  zero). 
In  general,  however,  a  detonation  implies  decrease  of  specific 
volume  (or  increase  of  density),  while  combustion  Implies  de- 
crease of  pressure. 

We  assume  the  nature  of  the  gases  is  such  that  p.  >  p  ; 
in  other  words,  that  a  constant  volume  detonation  raises  the 


»  The  proof  is  implied  in  relation  (4), 


CONFIDENTIAL 


C0UPID3NTIAL 


pressure  (cf«  Becker  [6]).   As  a  consequence  we  have  also 
fg  >  r  ;  In  other  words,  a  constant  pressure  deflagration 
raises  the  specific  volume.   We  further  assume  that  dp/dr<  0, 
d^p/dT  ^  >  0. 

To  exhibit  how  the  gas  velocities  behave  behind  a  reaction 
front  we  represent  the  transitions  by  a  graph,  representing  for 
fixed  u^,  p^,  r^  in  a  (u,p) -plane  the  relation  between  p  and 
the  velocity  jump  u  -  u^^  =  v  -  v^.   This  relation  is  obtained 
by  elimination  of  r    from  relations  IV  and  (6).   The  graph  in 
the  (u,p) -plane  is  indicated  in  Fig.  2. 


Figure  2 

Relation  between  pressure  and  velocity 

behind  a  reaction  front. 


d£ 
d- 


PV^ 


condition 


*Since  by  IV,  with /T  =  (^-iytJ?+l) 
dp/dr  <  0  is  equivalent  withr>/^2'r^  throughout,  which  we  shall 

assume.   Note  that  for  most  gases  "i"  -*■  1  and  hehce  /^^ 
as  9  ->  oo  . 
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A   transition  to  a   state  represented  by  a  point  on  the 
upper  branch  corresponds   to  a  detonation,  while  the  lower 
branch  corresponds   to  a  deflagration.      If  the  reaction  front 
Is  "facing  to  the  right,"  l»e.    If  the  explosive  gas    (0)   Is 
to   the  right. and  the  burnt  gas   to  the  left,    (the  positive  x-axls 
defining  the  direction  to  the  right),   the  mass  flux  Is   from  the 
right  to  the   left,   v     <  0,   and     consequently 


P-P, 


_o 

u-u. 


Therefore,  the  right  upper  and  the  left  lower,  branch  correspond 
to  forward  facing  reactions.   It  Is  then  clear  that  for  a  for- 
ward facing  detonation  u  >  u^,  while  for  a  forward  facing 
deflagration  u  <  u^.  In  other  words: 

A  detonation  accelerates  the  burnt  gas  toward  the  explosive, 
a  deflagration  accelerates  the  bxirnt  gas  away  from  the  explosive. 

Among  all  conceivable  reaction  processes  which  lead  from 
a  fixed  state  (p^,  u^,  r^)  to  states  (p,  u,  r,),  there  are  two 
of  particular  significance,  one  a  deflagration  and  one  a  detona- 
tlon.   They  are  characterized  by  the  " Chapman- Jouguet  law 
which  can  be  stated  in  various  equivalent  forms;  for  most  re- 
action processes  it  becomes  plausible  by  many  considerations  (see 
e.g.  subsequent  discussions  in  this  supplement),  but  must  be  con- 
sidered as  a  hypothesis  rather  than  a  law  of  general  validity. 
The  "Chapman- Jouguet  reactions"  can  be  defined  by  any  of  the 
following  conditions  P,  Q,  R,  S,  T: 

(P)  ^  =  ^^ 

^^'  .dr   r-r 


(That  there  are  Just  two  solutions  of  this  equation,  one  with 
p  >  p  and  one  wlthT>r  ,  follows  from  the  assumptions 
d^p/d-r^  >  0,  ^  <  0,  and  p^  <  p^  or  r^  <  Tg,  see  p.  10. )  On 
the  graph  of  p  as  a  function  of  r  (Pig.  1)  the  point  corresponding 

»  See  Manual  [1]. 
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to  this  transition  is  characterized  by  the  condition  that 
the  tangent  at  this  point  passes  through  the  point  (t  ,p  ). 
Similarly,  on  the  graph  In  the  (u,p) -plane  (Fig. 2)  this 
point  Is  characterized  by  the  condition  that  the  tangent 
at  this  point  passes  through  (u^,p^)*  .   Equivalent  properties 
are:**  (See  Manual  p.  131) 

(Q)  The  propagation  speed  |v^|  of  the  front  observed  from  the 
explosive  is  an  extremum  for  fixed  values  of  ^  fPojT  )  (Chapman 
1899),  or 

dv  =  0  . 

(R)   The  entropy  of  the  burnt  gas  Is  an  extremum. 

(S)  The  propagation  speed  Iv|  of  the  front  observed  from  the 
burnt  gas  Is  sonic  (Jouguet  1905) 


(T)  ^  =  -i'fi 

dr     r 


The  two  Chap man- Jouguet  reaction  processes,  corresponding 
to  the  points  D  and  C,  differ  In  that  for  detonations  the  speed 
(▼qI  and  the  entropy  are  relative  minima,  while  for  deflagrations 
the  speed  |v^|  and  the  entropy  are  relative  maxima,  the  state 
(0)  ahead  of  the  wave  always  assTimed  to  be  fixed  (see  App.I^).*** 

♦   Prom  (4)  we  have  by  (P) 

2(u-u^)du  =  ('r^-r)dp-cp-pQ)dT  =  2(r^-T)dp  ,  whence 
du/dp  =  (u-u^)/(p-p^). 

♦♦  (For  the  equivalence  of  the  conditions  (P)  to  (T)  (see  App.I.]0< 
A  procedure  for  computing  the  Chapman -Jouguet  reaction  process 
when  the  state  In  front  of  it  is  given  is  outlined  In  Appendix 

♦**  Emphasis  is  laid  on  the  relative  character  of  these  oxtrema. 
The  absolute  value  of  the  maximum  for  deflagrations  may  be 
less  than  that  of  the  minimum  for  detonations. 
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the  extremal  properties  of  the  speed  |v  |  can  be  tmderstood* 


from  the  relation  (4), 

2      2  _  P-Po 

P-Po 
and  the  fact  that      Is  the  slope  of  the  segment  connecting 

the  points  ('^q,Pq)  anS  ('r',p)  together  with  d  p/dt^  >  0. 

According  to  whether  or  not  a  discontinuity  effects  a 
smaller  or  larger  change,  a  detonation  implying  a  transition 
of  Pq  to  p  <  p  or  p  >  p_j  will  be  called  weak  or  strong 
respectively;  a  deflagration  implying  a  transition  of  p  to  p  >  p« 
or  p  <  Pj^  will  also  be  called  weak  or  strong  respectively.  Both 
a  constant  volvime  detonation  and  a  constant  pressure  deflagration 
are  weak. 

The  answer  to  our  main  question,  the  problem  of  deter- 
minacy  of  reaction  fronts,  will  depend  decisively  on  whether 
the  flow,  observed  from  the  reaction  front,  is  sonic  or  sub- 
sonic or  supersonic.   This  question  is  answered  by  a  statement 
(Jouguet  [5] (which  we  shall  prove  in  Appendix  1.3),  and  which 
will  serve  as  the  basis  of  the  subsequent  analysis: 
The  gas  flow  relative  to  the  reaction  front  is 

supersonic  ahead  of  a  detonation  front, 
supersonic  behind  a  weak  detonation  front, 
subsonic  behind  a  strong  detonation  front, 
subsonic  ahead  of  a  deflagration  front, 
subsonic  behind  a  weak  deflagration  front, 
supersonic  behind  a  strong  deflagration  front. 
We  shall  refer  to  and  make  use  of  these  properties  as  Jouguet 's 
rule  in  the  following  discussion. 

*  One  might  also  use  the  relation  (7) 

V  -      P'Pq 
o  =  "  '^o  u-u^,  ' 

p-p 
and  the  fact  that       is  the  slope  of  the  segment  connacting 

the  points  (u^jp^)  anS  (u,p).   This  fact  leads  to  an  obvious 

geometric  construction  of  v  in  the  (u,p)  -plane,  see  Courant- 
Priedrichs  [8],  p.  51.     ° 

COHPIDENTIAL 


11  CONFIDENTIAL 


2.  general  remarks  on  the  determinacy  of  a  flow  Involving 
a  discontinuity  front. 

As  was  said  before,  we  restrict  our  discussion  to  the 
case  of  gases  filling  a  long  tube  and  asstime  that  all  processes 
are  one -dimensional,  (an  assumption  which  will  be  actually 
satisfied  for  combustion  processes  only  If  the  tube  has  a  con- 
siderable width). 

Assuming  that  the  state  of  the  gas  at  some  Initial  time 
is  given  and  that  at  this  initial  time  a  reaction  front  occxxrs 
at  some  place,  we  ask  how  far  the  flow  is  then  determined  by 
the  three  laws  of  conservation  and  the  boundary  conditions  Im- 
posed at  the  ends  of  the  tube.   (These  conservation  laws  are, 
as  said  before,  nothing  but  the  flow  differential  equations  In 
the  continuous  regions  of  the  flow  and  the  transition  relations 
I,  II,  III  at  a  reaction  front. ) 

We  specify  our  initial  conditions  by  assvuning  that  at 
time  t  =  0,  the  explosive  is  at  rest  in  the  half  infinite  tube 
X  ^  0,  that  a  reaction  begins  at  t  =  0  at  the  closed  end  x  =  0, 
and  that  the  closed  end  is  operated  as  a  piston  moving  with  the 
prescribed  velocity  TJ. 

Usually  one  would  be  concerned  with  the  case  U  =  0,  which 
means  that  the  end  of  the  tube  is  held  fixed;  however,  consider- 
ing the  general  case  of  an  arbitrary  value  TJ  not  only  contributes 
much  to  the  clarification  of  the  problem  of  determinacy,  but  al- 
so is  helpful  for  the  understanding  of  phenomena  such  as  occur 
when  an  explosive  gas  moves  against  a  wall  and  is  ignited  there; 
likewise  for  an  understanding  of  the  phenomena  of  boosting,  the 
study  of  the  case  U  ^  0  may  be  useful. 

The  mathematical  determinacy  of  flows  involving  reaction 
fronts  depends  essentially  on  the  geometrical  relationship  be- 
tween the  characteristics  of  the  system  of  differential  equations 
I,  II,  III,  and  the  line  W  representing  the  reaction  front  in 
the  (x,t)  plane  as  well  as  the  line  ir  representing  the  moving 
piston  at  the  closed  end  of  the  tube.   As  usual,  the  one  dimensional 
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Figure  3 
Possible  relation  between  characteristic 
and  reaction  front*   (For  other 
possibilities  see  Flg.7.) 

motion  is  represented  in  the  (x,t)  -plane.  For  Isentroplc  flow 

the  system  reduces  to  the  system  of  the  differential  equation 

I  and  II  only  (see  Manual  [l].  Chapter  III)  and  for  any  flow 

there  exists  two  sets  of  characteristics  C_,  C^  In  the  (x,t)-plane, 

defined  (see  Manual,  Chapter  IV  )  by 

C^  :  dx  =  (u+c)dt;  along  C_|_  we  have  dp  +/?c  du  =  0 
C_  :  dx  =  (u-c)dt;  along  C_  we  have  dp  -/3c  du  =  0 

If  the  flow  Is  not  Isentroplc  the  third  family  of  characteristics 
C  ,  the  streamlines,  have  to  be  considered 


dx 


udt;  along  C  we  have  d  E  +  pdt 


For  the  following  discussion  we  envisage  the  case  of  isentroplc 
flow,  keeping  in  mind  that  the  situation  is  analogous  In  the 
general  case.  Furthermore,  if  nothing  is  said  to  the  contrary, 
we  shall  visualize  the  flow  in  front  and  behind  the  reaction 
front  as  consisting  of  zones  of  constant  flow  and  siniple  waves, 
so  that  either  both  sets  C_  and  C^  or  at  least  one  is  straight. 
(The  following  statements,  however,  remain  valid  also  for  flows 
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In  which  the  characteristics  are  all  curved).   We  Imagine 
the  characteristics  described  In  the  direction  of  Increasing 
values  of  t.   Then  a  line  element  X    through  a  point  P  In  the 
(x,t)  -plane  Is  called  space-like  If  the  two  characteristics 
leading  Into  P  are  both  on  the  same  side  of  ?i  ;      otherwise 
Is  called  time-like  (see  Fig.  3). 


Flgtire  4 


space-like 


time -like 


A  line  or  curve  through  a  zone  of  flow  In  the  {x,t)  -  plane  Is 
called  space-like  if  all  Its  line  elements  are  space-like. 
(E.g.  t  =  0,  the  x-axls  Is  space-like.) 

If  Initial  values  of  u  and  /O  are  arbitrarily  prescribed 
along  a  space-like  curve  L  (see  Pig.  5),  then  the  initial  value 


M-f 


Figure  5 
Domain  of  dependence, 
AB,  of  a  point  P. 
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probleia  Is  uniquely  solvable,  and  the  values  of  u  and  p   at 
a  point  P,  as  In  Pig.  5,  depend  only  on  the  Initial  values  of 
u  and  (o   on  that  part  AB  of  L  which  Is  cut  out  hy  the  two 
characteristics  drawn  backwards  from  L  domain  of  dependence). 

If  L  were  time -like,  then  only  one  of  the  two  char a ot er- 
istics would  Intersect  L;  along  a  time -like  line  not  both 
quantities  u  and  ^   can  be  prescribed  arbitrarily,  but  only 
one  of  them.  The  same  Is  true  If  L  Is  characteristic. 

If  u  and  (^    are  prescribed  on  the  space-like  Initial 
curve  L:  AB  then  the  flow  Is  determined  between  L  and  the 
characteristic  C°  through  A.  To  determine  the  flow  beyond 
0°  further  conditions  must  be  given.  In  particular,  let  us 
consider  the  sector  G  between  C^  and  a  time-like  curve  J,  then 
the  existence  of  one  solution  ^s  uniquely  determined  If  along 
G  only  one  of  the  quantities  u  or  /o  is  (arbltrarllij)  prescribed. 


Figure  6 
Along  space-like  L  two  quantities 
u,x7  can  be  prescribed;  along  time- 
like J  only  one  quantity. 
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Generally  speaking:   The  nvmb^r  of  conditions  to  be 
prescribed  along  a  portion  of  the  boundary  of  the  domain  In 
the  (x,t)  -  plane  Is  given  by  the  nxmber  of  Intersections 
with  characteristics  drawn  backwards  from  points  near  the 
considered  portion  of  the  boundary. 

In  the  case  of  a  non-lsentroplc  flow,  also  the 
third  set  of  characteristics  Cq,  the  streamlines, 
has  to  be  considered.  But  the  situation  remains 
similar.   We  draw  from  a  point  P  the  characteristics 
backwards  (towards  decreasing  values  of  t ) j   then 
along  a  space-like  part  of  the  boxmdary  of  a  region 
g  in  the  (x,t)  -  plane  three  unknown  quantities,  u, 
fo    and  say  the  pressxire  p,  can  be  prescribed.   The 
munber  of  quantities  prescribable  along  a  time-like 
line  T  is  one  if  only  one  of  the  characteristics, 
say  C+,  through  points  of  J  leads  Into  g,  while 
two  quantities  may  be  prescribed,  when  there  are 
two  such  characteristics,  say  C_^  and  C  . 

After  these  remarks  the  various  possibilities  of  flows 
Involving  reaction  fronts  can  easily  be  analyzed.   The  domains 
G  in  which  we  are  primarily  interested  are  first  that  between 
the  (space -like)  x-axis  and  the  reaction  line  W  representing 
the  front  in  the  (x,t)-  plane,  and  second  the  region  between 
the  reaction  line  W  and  the  piston  curve  tt  :  x  =  Ut. 

We  distinguish  four  cases  AA,  AB,  BA,  BB  according  to 
whether  the  flow  relative  to  the  line  W  is  supersonic  (A)  or 
subsonic  (B)  before  or  behind  the  front.   In  a  schematic  way, 
visualizing  W  and  the  characteristics  as  straight,  these  cases 
are  represented  by  the  following  diagrams: 
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Case  AA 


Case  A5 


Case  BA 


Case  BB 


Figure  7 
Possible  relations  between  characteristics 
and  reaction  front. 
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TMs  follov;e  frcm  the  expressions  for  the  directions  -^  along 
the  characteristics.   In  the  case  of  supersonic  flow  in  front 
of  the  wave  line,  W  is  space-like  observed  from  the  region 
ahead  of  W;  hence  the  initial  values  along  the  x-axls  at  t  =  0 
determine  the  flow  uniquely  up  to  W,  no  matter  what  the  speed 
of  the  wave  front  is,  as  long  as  the  flow  ahead  of  W  remains 
supersonic  relative  to  W,   Then  by  the  transition  conditions  IV, 
V,  VI,  the  values  of  p,  u,  /?   immediately  behind  the  wave  front 
are  determined.   Now  there  are,  for  any  such  wave  velocity,  two 
cases.   In  the  first  case,  AA,  also  the  flow  behind  W  is  super- 
sonic relative  to  W.   Then  W  is  space-like  for  the  region  behind 
W,  and  in  the  domain  between  the  wave  front  W  and  the  piston 
curve  the  flow  is  uniquely  determined  by  the  additional  boundary 
condition  u  =  U  at  the  piston  (assuming  the  piston  curve  time- 
like). 

In  the  second  case,  AB,  again  the  state  ahead  of  W  is 
uniquely  determined  by  the  initial  conditions  along  x  ^  0  for 
t  =  0.   Then  the  transition  conditions  determine  the  state  Imme- 
diately behind  W.   How,  however,  W  Is  not  space-like  with  respect 
to  the  zone  behind  W:   Therefore  the  state  immediately  behind  W, 
considered  as  initial  state  for  the  region  between  W  and  it,    can  no 
longer  be  prescribed,  but  has  to  be  subjected  to  a  restrictive 
condition  to  become  compatible  with  the  bo\mdary  condition  at  the 
piston  curve  tr.   While  thus  in  the  case  AA  the  wave  front  W  could 
be  chosen  arbitrarily,  there  is  no  leeway  In  the  choice  of  W  in 
the  case  AB.   The  situation  is  similar  to  that  for  non-reactive 
flow  involving  a  shock  line  W;   the  velocity  of  the  discontinuity 
front  is  determined  in  such  a  way  that  adjustment  to  the  piston 
motion  results. 

In  the  case  BA  the  wave  line  W  observed  from  the  state 
ahead  of  it  is  time-like,  so  that  in  front  of  W  a  boundary  condi- 
tion can  be  arbitrarily  imposed.   Consequently  the  wave  speed 
and  one  quantity  ahead  of  the  wave  could  be  prescribed:  in  other 
words,  there  are  two  degrees  of  indeterminacy  in  the  case  BA. 
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Thus,  In  addition  to  the  Initial  state  and  the  piston  motion 
two  more  quantities  can  and  must  "be  prescribed  to  determine 
the  flow. 

In  case  BB,  only  one  quantity  such  as  the  wave  velocity 
can  be  prescribed;  the  other  quantity  In  front  of  the  wave 
must  then  be  so  chosen  that  the  resulting  values  of  the  three 
quantities  behind  the  wave  are  compatible  with  the  boundary 
condition  at  the  piston. 

It  should  be  emphasized  that  in  both  cases  BA  and  BB 
the  state  ahead  of  the  wave  is  Influenced  by  the  wave* 

3«  One-dlmenslonal  gas  flow  involving  a  reaction  front. 

Combined  with  Jouguet's  rule,  the  preceding  remarks  now 
permit  a  complete  discussion  of  our  problem  of  determlnacy  for 
gas  flow  Involving  a  detonation  or  deflagration  front.   The 
higher  degree  of  indeterminacy  of  such  flows  as  compared  with 
flows  involving  merely  shock  discontinuities  is  due  to  the  fact 
that  all  cases  AA,  AB,  BA,  BB  may  occur  in  reaction  processes, 
while  shocks  always  belong  to  the  case  AB.   In  the  following 
detailed  discussion  we  desist  from  a  systematic  eniimeration  of 
all  mathematical  possibilities  of  flows  compatible  with  the 
conditions,  and  instead  consider  the  most  typical  and  important 
cases* 

Since  the  state  (0)  ahead  of  the  front  is  considered  as 
given,  the  reaction  wave  i.e.  Its  speed  and  the  state  (l)  imme- 
diately behind  it,  is  determined  by  one  parameter  for  which  we 
choose  the  velocity  u  behind  the  front.   The  speed  of  the  wave 
then  is  given  by  v^  =  -~  2^a  .  we  first  assume  that  W  represents 
a  detonation,  with  a  state  of  rest,  u^  =  o,  ahead  of  it.   Then 
as  seen  before  we  have  u  >  0.  By  Uj.  =  Cj.  we  denote  the  parti- 
cular value  for  the  gas  velocity  behind  the  detonation  front 
which  satisfies  the  Chapman- JoTiguet  condition. 

We  start  with  the  discussion  of  a  weak  detonation,  which 
means  u  <  Uj^;  according  to  Jouguet's  rule  this  corresponds  to 
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the  cases  AA,  and  W  Is  space-like  with  respect  to  both  adjacent 
regions.  As  a  consequence.  In  view  of  the  statements  In  the 
preceding  section,  every  weak  detonation,  i.e.,  every  speed  for 
the  detonation  front  compatible  with  the  condition  u  <  Uq  Is 
mathematically  con5)atlble  with  the  data.   The  state  of  the  flow 
In  front  of  W  is  determined  by  the  Initial  condition,  the  state 
immediately  behind  W  by  the  transition  conditions,  and  It  remains 
to  analyze  In  detail  how  the  adjustment  of  the  velocity  of  the 
flow  behind  the  front  to  the  piston  velocity  Is  achieved  for  a 
fixed  choice  of  our  parameter  u.   This  adjustment  depends  on 
whether  or  not  the  piston  speed  U  exceeds  the  value  Uj^: 

Case  1.   U  <  Uq.   In  this  case  every  weak  detonation,  u  chosen 


arbitrarily  in  the  interval  0 


is  possible. 


We  first  consider  detonations  with  U  »  u  «  u-^.   Then  the 
transition  from  the  flow  with  u  >  U  behind  the  front  to  a  flow 
with  u  =  U  at  the  piston  is  effected  by  a  centered  rarefaction 
wave  (see  Pig.  8). 


/P<a  re  fa C  tio  ?7 


Figure   8 

Weak  detonation  with  U  <  rx  <  u-^ 

for  piston  with  U  <  u-^. 
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In  the  limiting  Chapman- Jouguet  case  u  =  Uj>,  which 
separates  weak  and  strong  detonations,  the  rarefaction  wave 
follows  the  detonation  wave  immediately  (see  Pig.  9): 


Ro.ref -Action 


onat/or? 


Figure  9 
Detonation  with  U  <  u  =  Uj.  satisfying 
Chapraan-Jougiiet»s  hypothesis. 


In  case  u  =  U,    the   rarefaction  wave   drops   out* 
If   the  piston  speed  u  exceeds  u: 

0  <  u  <  U 
the  adjustment  is   effected  by  a   shock    (see  Pig.    10). 

CONFIDENTIAL 


21 


CONFIDENTIAL 


Plgtire  10 
Detonation  with  u  <  TJ  for  piston 
with  u  <  Mj^» 


Summing  up:   If  U  <  u^  a  variety  of  weak  detonations  Is  compatible 
with  the  condition  of  the  problem;  specifically,  any  value  of 
the  velocity  u  behind  the  detonation  front  ranging  from  u  =  0 
(constant  volume  detonation  with  Infinite  detonation  speed)  to  the 
value  u  =  Uj.  ( Chapman- Jouguet)  Is  possible. 

It  should  be  mentioned  that  there  Is  an  additional  mathemati- 
cal Indeterminacy  In  the  following  sense:  At  any  time  after 
the  start  of  the  process  the  wave  may  change  Its  character; 
then  It  may  continue  as  a  stronger  (though  still  weak)  detona- 
tion, followed  by  a  rarefaction  wave  which  will  have  to  be 
balanced  by  a  backward  rarefaction  wave,*  or  It  continues  as 
a  weaker  detonation  wave,  followed  by  a  shock  which  Is  balanced 
by  a  backward  shock  (see  Pigs. 11, 12). 


♦  As  can  be  seen  by  the  method  described  In  the  Manual  [1], Chap. Ill, 
Section  48. 
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Figure  11 

Detonation  changing  to 

a  weaker  one* 


Figure  12 
Detonation  changing  to 
a  stronger  one. 


Likewise,  as  appears  from  these  last  remarks,  non- 
constant  gradually  changing  detonation  waves  are 
compatible  with  the  conditions* 


Case  2.   U  >  Uj..   Also  in  this  ease  a  variety  of  weak  detona- 
tions are  possible;   the  adjustment  to  the  piston  velocity 
Is  then  always  effected  by  a  shock.  Again  the  range  of  these 
detonations  can  best  be  described  by  varying  the  parameter 
u  from  u  =  0  (constant  volume  detonation)  through  positive 
values  and  by  using  a  (p,u) -diagram. 

It  can  be  seen  that  there  is  a  value  u  =  u„  of  the  para- 

T   p"p_ 
meter  with  u^  <  u^j  for  which  the  velocity  ^ ^  ^  ±-       ^     of  the 

weak  detonation  is  eqtial  to  the  velocity  of  the  shock  which 
accelerates  the  gas  to  the  piston  speed  U  >  u^..   How  every  weak 
detonation  with  u  <  u  is  possible  followed  by  a  (slower)  shock 
which  adjusts  the  particle  velocity  to  the  piston  speed  U.   If 
the  particle  velocity  u  behind  the  weak  detonation  front  Is  equal 
to  u  then  the  detonation  front  and  the  subsequent  shock  coalesce 
Into  one  single  discontinuity  front  W.  Across  this  combined 


COHFIDENTIAL 


29 


COHFIDEHTIAL 


front  the  laws  of  conservation  are  tatlsfled  and  the  pressure 
la  ralsedj  hence  the  combined  front,  represented  by  W,  Is  a 
detonation  front*   Since  the  velocity  behind  the  combined 
front  l3  u  =  U  >  rx~,   this  detonation  Is  a  strong  detonation 
front,  appearing  as  a  weak  detonation  Immediately  followed  by 
a  shock* 

While  In  the  case  TJ  >  «_.  there  are  again  many  weaic  de- 
tonations (with  u  <  u  )  possible,  there  exists  one  well  deter- 
mined strong  detonation,  directly  producing  a  constant  flow  of 
the  same  speeds  as  the  piston  (see  Pigs ,13  and  14).   This  Is 
In  agreement  with  the  fact  that,  according  to  Jouguet*s  state- 
ment, a  strong  detonation  cerresponds  to  the  case  AB* 


Flgtire  13 
Detonation  for  piston 
with  u  >  u_. 


Figure  14 

Detonation,    for  piston  with 

u  >  u_.    Immediately  followed 

by  a  shock, both  coalescing 

Into  a   strong  detonation. 
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As  In  Case  I,  the  continuation  of  the  wave  is  not 
uniquely  determined;   later  changes,  as  explained 
In  that  case,  are  coatpatlble  with  the  conditions 
of  the  problem. 

A  strong  detonation  for  XJ  <  Uj.  is  not  compatible  with 
the  conditions;   again  we  would  have  the  case  AB;   for  the 
domain  behind  the  line  W  both  W  and  the  piston  line  ir  are  time- 
lifce  and  along  each  of  these  lines  u  is  constant;  but  no  solu- 
tion of  the  differential  equations  between  two  time-like  lines 
can  have  such  boundary  values  except  the  solution  u  «  const; 
since  U  <  Ujj  at  tf  and  u  >  Uj^  at  W  we  have  a  contradiction,  which 
shows  that  no  adjustment  of  a  strong  constant  detonation  to  a 
slow  piston,  e.g.  a  fixed  wall,  is  possible. 

As  we  have  seen,  however,  a  Chapman- Joup;uet  detonation 
can  always  be  adjusted  to  the  piston  motion. 

4,  Remarks  on  the  Chapman- Jouguet  hypothesis. 

Taking  into  account  the  physical  situation,  Jouguet  has 
advanced  the  following  arg\ment  in  support  of  the  Chapman- Jouguet 
condition. 

1«   If  the  detonation  is  to  proceed  with  a  certain  velocity, 
there  are  still  two  detonations  possible,  a  weak  one  and  a  strong 
one,  the  latter  equivalent  to  the  weak  detonation  followed  by  a 
shock.   The  strong  detonation  leads,  therefore,  to  the  higher  en-  - 
tropy  of  the  burnt  gas.  Since  the  state  with  higher  entropy  is 
the  more  likely  one,  most  probably  the  strong  detonation  will  occur. 

2.  Due  to  heat  conduction  and  rupture  of  the  enclosure  of 
the  explosive,  a  rarefaction  wave  will  be  sent  after  the  detonation 
wave. 

3.  Since  the  flow  behind  the  strong  detonation  front  is 
subsonic  relative  to  the  front,  the  rarefaction  wave  will  catch 
up  with  the  detonation  and  weaken  it  until  the  flow  behind  the 
front  has  become  sonic,  in  other  words,  until  the  Chapman- Jouget 
condition  is  satisfied. 

♦  See  Becker  [6],  p.  354.   , 
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Objection  might  be  raised  to  the  first  step  of  the 
argument.  Whether  a  weak  detonation  or  a  weak  detonation 
immediately  followed  by  a  shock  (I.e.  a  strong  detonation) 
occurs,  should  depend  on  whether  or  not  the  resulting  flow 
la  compatlblo  with  the  conditions.    These  conditions  are 
not  explicitly  stated  by  Jouguet  and  Becker,  but  a  fixed 
closed  end  seems  to  be  assumed.   Then  as  seen  before  a  strong 
detonation  Is  not  compatible  with  the  conditions. 

However,  one  may  argue  Instead  that  the  start  of  a  de- 
tonation is  a  complicated  process  which  cannot  properly  be 
described  by  gas  flow  involving  a  single  discontinuity  front. 
If  such  a  description  is  appropriated  for  a  somewhat  later 
instant,  adjustment  of  a  s trong  detonation  wave  to  a  closed 
end  could  be  achieved  by  a  rarefaction  wave  sent  back  toward 
the  closed  end  and  being  reflected  there.   This  reflected  wave 
would  then  overtake  the  detonation  wave  and  weaken  it  until 
the  Chapman- Jouguet  condition  is  satisfied.   (See  Step  3  of 
Jouguet ' s  argument )  • 

Thus  the  second  step  of  Jouguet* s  argument  appears  to 
be  unnecessary.  Although  rarefaction  waves  due  to  rupture 
and  heat  loss  will  occur,  a  rarefaction  wave  is  already  needed 
for  adjustment  of  the  flow  to  the  closed  fixed  end. 

The  final  theoretical  decision  about  the  Chapman- Jouguet 
hypothesis  is,  of  course,  in^ossible  without  detailed  investi- 
gation of  the  reaction  process  involved.   Experimentally  the 
Chapman -Jouguet  hypothesis  seems  to  be  well  confirmed.   There- 
fore it  appears  reasonable  to  accept  the  hypothesis  for  calcula- 
tions, even  when  non-steady  processes  are  involved.   Numerical 
procedure  is  facilitated  because,  as  we  saw,  the  flow  behind 
a  constant  Chapman- Jouguet  detonation  is  a  centered  rarefaction 

♦  Just  as  the  decision  as  to  whether  or  not  a  shock  occurs 
depends  on  whether  or  not  the  resulting  flow  is  compatible 
with  the  conditions  of  the  problem  and  not  on  the  fact  that 
the  entropy  increases  through  a  shock. 
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wave,  the  first  ebaracterlstlo  of  which  oolncldea  with  the 
detonation  front*   This  fact,  apparently  not  mentioned  In 
the  literature  explicitly  but  Infilled  In  Taylor » a  work  (see 
[14]),  Is  decisive  for  actual  calculation. 

If  a  detonation  process  Is  to  he  calculated  for 
which  the  medium  ahead  of  the  wave  Is  not  In  a 
constant  state,  and  If  the  Chapman-Jouguet  hypoth- 
esis Is  assumed  to  be  valid,  then  one  should  be 
aware  of  the  fact  that  the  detonation  front  Is  no 
longer  a  characteristic  of  the  flow  behind  It;   It 
Is  rather  an  envelope  of  (not  necessarily  straight) 
characteristics  (see  Pigs,  15,16). 


Figure  15 
Accelerating  Chapman- Joviguet 
detonation  wave* 


Figure  16 
Decelerating  Chapman- Jouguet 
detonation  wave. 
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(In  specific  oases  It  was  found  by  conrputatlon 
that  the  deviation  of  the  detonation  front 
from  a  single  characteristic  Is  very  slight.) 

5,   In  deflagration  processes  the  situation  Is  In  many 
respects  quite  different  from  that  encountered  with  detona- 
tion processes.   Suppose  a  weak  deflagration  wave  (practi- 
cally all  observed  deflagration  waves  seem  to  be  weak)  be- 
gin at  the  piston,  x  =  0,  and  moves  Into  the  unburnt  gas  In 
the  interior  of  the  tube,  x  =  0»     Then  the  velocity  u  of  the 
burnt  gas  behind  the  deflagration  front  Is  negative,  u  <  0. 
This  is  compatible  with  the  conditions  of  the  problem  only 
if  the  piston  Is  moved  out  with  a  speed  at  least  equal  to  u. 
Otherwise  It  Is  In^josslble  that  a  deflagration  front  moves  In- 
to the  explosive  gas.  For,  we  have  here  the  case  BB;   the 
line  W  is  time-like  with  respect  to  the  region  behind  it;  an 
adjustment  to  a  forced  or  forward  moving  piston  is  impossible. 
What  one  observes  actually,  is  that  a  pre-compression  wave  is 
sent  out  Into  the  explosive*   It  pushes  the  explosive  ahead 
with  a  Telocity  just  sufficient  to  insure  that  it  may  come  to 
rest  when  it  is  swept  over  by  the  deflagration  wave* 

The  occurrence  of  a  pre-compresslon  wave  is  in  complete 
agreement  with,  or  rather,  a  consequence  of,  the  theoretical 
considerations  of  Section  2.   Indeed,  according  to  Jouguet'a 
rule,  the  flow  ahead  of  a  deflagration  is  subsonic  (case  B) 
and  consequently  the  deflagration  Influences  the  state  of  the 
gas  ahead  of  it.   Observation*  shows  that  the  deflagration 
wave  moves  with  gradually  increasing  velocity  and  sends  out 
a  continuous  coii?)resslon  wave  into  the  explosive.  A  sin^jlified 
model,  however,  is  more  convenient  for  theoretical  treatment: 
to  assume  that  a  single  constant  shock  wave  is  sent  out  into 
the  explosive  and  that  the  deflagration  wave  follows  with  a 
constant  velocity  (see  Pigs.  17  and  18). 

*  See  Payman  [9],  [ lOj ,  [11],  [12],  [13]. 
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Figure  17 
Actual  deflagration  with 
continuous  pre-conipresslon, 
U  «  0. 
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Figure  18 
Deflagration  with 
shock  pre-compresslon, 
U  »  0. 


To  each  assumed  pre-con5)reaslon  shock  there  is  still 
a  great  variety  of  possible  adjustments  either  by  a  uniquely 
determined  weak  deflagration  wave  (case  BB)  or  by  a  strong 
deflagration  wave  (case  BA)  which  may  be  arbitrarily  chosen 
within  a  one  parameter  variety,  and  Is  followed  again  by  rare- 
faction or  shock  waves  (see  FlgSf  19  and  20)* 
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Figure  19 
Deflagration  followed 
by  shook. 


Figure  20 
Deflagration  followed 
by  rarefaction. 


These  possibilities  are,  of  coxorse,  limited  by  the  condi- 
tion that  the  Telocity  of  each  wave  is  less  than  that  of 
the  preceding  wave.   If  the  shock  wave  following  a  strong 
deflagration  wave  coalesces  with  it  the  result  is  equivalent 
to  a  weak  deflagration  wave.   A  deflagration  wave,  however, 
coalescing  with  its  pre-coBg>resslon  wave  is  equivalent  to  a 
detonation  wave.   Thus,  by  Increasing  the  strength  of  the 
pre -congress ion  wave  so  much  that  its  speed  equals  that  of 


•  Cf«  Jouguet  [5],  V.  Veumann  [7]. 
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a  possible  detonation,  one  may  envision  a  continuous  transi- 
tion from  deflagration  processes  to  detonation  processes. 

The  main  point  In  this  discussion  Is  that  deflagra- 
tlon  processes  have  a  higher  degree  of  Indeterminacy  than 
detonation  processes.  The  Chapman- Jouguet  condition  does  not 
make  the  process  unique;  for  sufflclfintly  strong  arbitrary 
pre-compresslon  shocks  a  rarefaction  wave  following  the  caiap- 
man- Jouguet  deflagration  can  always  be  found  such  that  the 
flow  at  the  piston  has  the  right  velocity.  To  determine  the 
process,  new  physical  conditions  will  have  to  be  added  that 
mast  arise  from  a  study  of  coiribustlon  processes. 

It  Is  clear  that  there  are  still  more  possibilities  in 
addition  to  those  mentioned,  viz.  those  in  which  the  pre-com- 
presslon Is  not  achieved  by  a  shock,  but  by  a  continuous  com- 
pression wave.  As  was  mentioned  before,  this  is  what  actually 
happens  when  deflagration  begins  at  one  end  of  the  tube.  The 
Chapman- Jouguet  condition  can  then  not  be  satisfied  in  the 
first  stage  of  the  process,  but  only  one  additional  condition 
Is  needed  as  long  as  the  deflagration  is  weak.  When  the  de- 
flagration becomes  strong  a  second  condition  (such  as  the 
Chapman- Jouguet  condition)  must  be  added. 

6*  Stationary  deflagration  fronts. 

So  far  we  have  mainly  been  concerned  with  reaction 
waves  moving  In  to  an  ejcplosive  gas  at  rest.  It  is  nat\irally 
of  great  Interest  to  study  other  situations  such  as  flows  in- 
volving a  stationary  reaction  front. 

Let  us  consider  In  particular  the  case  that  a  steady 
flow  of  contoustlble  gas  enters  one  end  of  a  cylindrical  tube, 
undergoes  combustion  on  passing  through  a  deflagration  or 
"flame"  front,  and  escapes  as  burnt  gas  at  the  other  end  of 
the  tube.   (In  a  crude  approximation  the  flow  through  the 
combustion  chasiber  of  a  jet  engine  might  be  described  in  such 
a  fashion. ) 
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Suppose  the  velocity  of  tbe  entering  explosive  gas 
is  varied,  the  state  of  the  explosive  being  kept  fixed; 
then,  if  at  all,  two  deflagration  processes  are  possible, 
a  weak  and  a  strong  one.  We  may  assume  the  weak  one  occurs. 
If  the  entrance  velocity  is  increased, the  pressure  drop 
across  the  deflagration  front  will  increase  until  for  a 
certain  value  of  the  entrance  speed  v-  the  Chapman- Jouguet 
reaction  is  reached.   A  f\u:»ther  increase  in  pressxire  drop 
would  make  the  deflagration  a  strong  one  and  at  the  same 
time  reduce  the  entrance  velocityj   for,  by  property  Q 
(see  p«  9),  the  relative  speed  between  explosive  and  reac- 
tion front  is  a  maximum  for  a  Chap man- Jouguet  deflagration. 
If,  therefore,  the  entrance  speed  is  increased  beyond  v  , 
a  stationary  deflagration  process  is  no  longer  possible, 
and  the  cojiibustion  will  therefore  be  stopped.   This  phenom- 
enon has  been  termed  "choking."   Its  occxirrence  is  another 
illustration  for  the  ingsortance  of  the  Chapman -Jouguet 
condition. 

7»   Remark  on  Taylor *s  Theory  of  Spherical  Detonation  Waves. 

Taylor  has  shown  that  the  flow  produced  by  a  spherical 
detonation  wave  is  accessible  to  mathematical  treatment.   The 
detonation  is  assumed  to  be  weak  or  to  be  a  Chapman- Jouguet 
detonation.   It  is  then  not  influenced  by  the  flow  behind 
it.   This  was  shown  earlier  to  be  true  for  one -dimensional 
waves  and  the  fact  remains  true  for  spherical  waves.  There- 
fore, it  is  consistent  to  assume  that  the  detonation  front 
proceeds  with  constant  radial  speed.   In  this  respect  a  spheri- 
cal detonation  wave  differs  from  a  spherical  blast  wave,  which 
is  influenced  by  the  flow  behind  it  and  gradually  fades  out 
unless  energy  is  constantly  added* 
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F«r  the  flow  behlod  the  wave  it  Is  assumed  that 
velocity  u,  pressiire  p,  and  specific  volxune  r  depend 
only  on  the  ratio  ^  =  p/t,  r  being  the  distance  from  the 
center  and  t  the  time*  This  as3uiiq)tion  is  consistent  with 
a  constant  speed  |  =  h:  of  the  detonation  front.   The  differ- 
ential equations  for  spherical  flow 

(8)  u^  +  u\iy  +  tpp  =  0,   r^  +  uTy  -  («^  •♦■  2r"^)r  «  0 

together  with  p"^dp  +  r'^dr  =  0  (see  Manual  [1],  V  90,p.263, 
(54),  (35)),  then  become 

(9)  (u  -  |)uw  -  '!?pr   =0     and 

(u  -0T-(u^+2r"^r=  0  , 

(of.  Manual  V90,  p.  265,  (34) »,  (35)'  )  .   Elimination  of 
Tv  and  introduction  of  c  =  "^pr  yields  the  equation 

(10)  u^  =  2k,'\x   [(i^)2  -  1]"^  . 

This  equation  was  solved  numerically  under  the  conditions  that 
behind  the  detonation  front  r  =  st,  the  quantities  u,  p,  r 
assvune  the  values  which  are  uniquely  determined  by  the  transi- 
tion conditions  of  the  reaction  process.   If  this  process 
satisfies  the  Chapman -Jouguet  condition,  the  rate  of  change  of 
velocity  with  respect  to  the  variable  ^  is  infinite  behind 
the  front  (while  it  Is  finite  for  one -dimensional  reactions). 
This  is  seen  from  (10)  since  condition  S  (p.  9  )  implies  that 
H  -  u  =  |v|  =  c«  If  the  detonation  is  weak,  the  rate  of  change 
of  velocity  is  finite  but  different  from  zero  In  contrast  to 
one -dimensional  flow  behind  a  weak  detonation  wave, which  im- 
plies a  section  on  which  the  state  is  constant  (see  Pig.  8). 
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However,  the  important  feature  of  tbe  solution  la  that  for 
a  certain  positive  value  of  r/t,  say  r/t  =  4-»  the  velocity 
Is  zero.   The  path  r  =  ^^t  Is  a  characteristic  which  separate 
the  rarefaction  wave  between  r/t  =  ^  and  r/t  =  H  from  a 
constant  state  between  r/t  =  0  and  r/t  =     K   • 

Thus,  a  spherical  flow  Involving  a  detonation  wave 
with  constant  speed  and  an  expanding  core  in  which  the  state 
Is  constant  and  the  gas  at  rest,  is  found  to  be  consistent 
with  the  general  laws  of  detonation. 

By  determining  such  a  flow,  Taylor  has  disproved  the 
statement  made  by  Jouguet,  quoted  on  occasion,  that  a  spheri- 
cal detonation  wave  could  not  maintain  itself. 
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Appendix  I 

In  this  Appendix  we  present  justifications  for  some 
statements  made  in  the  text. 

I,l.   We  first  prove  the  equivalence  of  the  various  forms 
(P)  to  (T)  of  the  Chapman- Jouguet  condition: 
Prom  (1)  and  (2)  we  have 

(V-  l)(dE  +  pdr)  =  T(dp  +'ir'\dT)   =  r(dp  -t-^^c^dr)} 
from  IV  and  (4)  we  have 

2(dE  +  pdr)  «=  (T^-r)dp  +  (p-PQ)dT 

«  (T^-T)  (dp  +  /oVdr) 


<^o-^)^/^H 


Vanishing  of  any  of  these  terms  is  therefore  equivalent. 
Relation  dE  +  pdr  «  o  expresses  that  the  entropy  is  station- 
ary (R)j  relation  dp/dr  =  -"^p/t  la  condition  (T)j  ^'sla- 

tlon  dp/dr  =  -(p-p„)/(r^-T)  is  condition  (P)}   relation 

2  o  '   o  2  2  2  2 

dv^  at  0  is  condition  Q.  Vanishing  of  dp  +/<rc  dr  and  dp+/0  v  dr 

is  equivalent  with  |v|  =  0  or  condition  (S). 

1.2.   To  establish  that  the  speed  |v  |  and  the  entropy  have 
a  minimum  aad  a  maximum  in  the  Chapman- Jouguet  states  D  and 
C,  we  should  show  that  the  derivatives  of  -^  and  of  9"-'-(^+p) 
with  respect  to  r  are  positive  in  state  D  and  negative  in 
state  C.  Since  these  expressions  vanish  at  these  points  and 
T  -  T  is  positive  in  D,  negative  in  C,  it  is  sufficient  to 
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show  that  the  derlratlve  of 

p-p 

Is  positive  In  D  and  C.   Since  M-:^r-^)   =  0  In  D  and  C,  the 
statement  follows  from  the  a3Sun?)tion  that  d  p/dr>  0  (see  p.10). 

1.3    In  order  to  Justify  Jouguet^s  rule  we  make  an  addition- 
al assximptlon  on  the  nature  of  the  burnt  and  unburnt  gases: 
Throuj^h  a  detonation  the  pres3\xre  rises  more  than  It  would 
through  adlabatlc  compression  to  the  same  volume,  and  through 
deflagration  the  pressvire  drops  less  than  through  adlabatlc 
expansion  to  the  same  voltune*  Let  p^,  t^,  p^^,  t  ^  be  pressure 
and  specific  volume  ahead  of  and  behind  the  reaction  front 
and  let  p  be  the  value  of  the  pressure  by  adlabatlc  change 
from  the  state  (p  ,r  )  to  a  state  with  T  »  r^.   Then  the 
assumption  above  Is  sln5)ly 


for  detonation  and  deflagration. 

We  further  employ  for  the  unbtirnt  gas  the  standard 
assun5)tlon  (see  Manual  I-2(a),  p.6) 

^2 

"S-X  >  0  for  adlabatlc  change 

d-r 

observing  that  the  sound  velocity  c^  Is  given  by 

we  have  as  a  consequence  of  the  latter  assumption,  depending 
,2  5  ^2  "^'^^ 


o         o  r  -r 
o  1 
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and  by  virtue  of  the  assumption  P^  '*^  Pi  » 


2  .  ^2  Pl"Po 


whence  by  (4), 


O   <  T 

o  >  o 


2  <  2 
o    o 


Since  T,  <  T  for  detonations,  ^^  *  '^^  ^^^   deflagra- 
tions we  have  proved  that  the  flow  ahead  of  the  reaction  front 
and  observed  from  it  is  supersonic  for  detonations  and  sub- 
sonic for  deflagration* 

To  prove  Jouguet^s  rule  for  the  state  behind  the  front 
we  first  obtain  from  the  relations  enumerated  in  App.  I.l  the 
expression 

7  =  ^*TiiFy^^-^'<^^*^^- 


The  derivative  with  respect  to  r  for  the  Chapman- Jouguet 
states  D  and  C  , 

according  to  the  assuniptlons  made  eturlier  (p»  10) •  The  Mach 
number  M  =  v/c  equals  1  In  the  states  D  and  C  by  condition  S, 
and  these  are  tiae  only  states  satisfying  the  Chapman- Jo\iguet 
condition  •  Therefore  M  >  1  if  'T'<  r  or  'r<  r  •   in  view 
of  dp/dr  <  0  this  is  equivalent  to  p  <  p^,  p  <  p^;  the  former 
case  corresponds  to  a  weak  detonation,  the  latter  to  a  strong 
deflagration  (see  p«  15 )•  Similarly  M  <  1  for  p  <  Pjj  and 
P  <PC  . 

*  This  argument  is  similar  to  that  for  shocks  in  the  Manual, 
see  III-32,p.  80.   If,  instead  of  the  assunption  p  <  p,, 
one  assumes  that  the  sound  velocity  increases  through 
a  constant  volume  detonation,  the  argument  of  H»  Weyl([l5] 
p.  10)  could  also  be  used* 
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Appendix  II 

Outline  of  a  procedure  to  eo]is?ute  a  CHiapman'-Jouguet 
reaction  If  the  state  In  front  of  It  la  given. 

We  shall  Indicate  briefly  how  to  determine  the 
Chapman- Jouguet  reactions  when  the  state  of  the  explosive 
is  given.   The  given  q\xantlties  are  pressure  p  ,  specific 
volujne  T  ,  the  quantity  9  ,  which  is  proportional  to  the 
temperature,  and  the  total  energy  per  unit  mass,  E  .  Further, 
the  function  E  =  E(©)  (and  consequently  V  =  V(©)  =  1  +  d9/dE) 
are  assumed  to  be  known.   The  problem  is  then  to  determine 
p,  ^  ,  ©,  E  on  the  basis  of  the  relations  I,  II,  IV  and  one 
of  the  conditions  (P)  to  (T). 

We  proceed  to  eliminate  v^,  v,  r  ,  and  /S  obtaining 
one  equation  for  the  quantity  0,  which  then  is  to  be  solved. 

Eliminating  v^  from  the  relations  I  and  II,  setting 
V  =  c  by  condition  S  and  c  =  V  ©  by  (2)  we  obtain 

^  "  Po  =  /"  <  -^  -  li  ^  ®  • 

Substituting  r  k  ©/p  and  r^  =  ®(/Po  ^^  *^  latter 
equation  and  In  IV  we  obtain 

P     © 

P     «  Po 
and 

E-E     -      ©  p 

Our  intention  Is  to  eliminate  P^/p  ;  before  doing 
this  we  Introduce  the  dimentionless  qtiantities 
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Then  the   two  equations  above  become 
2^^   =  1  -  t^ 

Eliminating    f    we  obtain 


^   =^     i(V+l)   \/(V+l)^^^-(^+l)^+2(l+y>/) 


2+^-2>7    -^^^ 


—  _. > 

1  (^+i)\/a+i)^-4^7?  +()? +1)^-2(1+ ^7?  ) 

the  upper  sign  corresponds  to  detonations,  the  lower  sign 
to  deflagrations. 

Suppose  9  and  E  are  given,  then  nq ,    V  ,  and  t,   depend 
on  9  only.  The  latter  equation  is  therefore  an  equation  in 
9  only,  which  can  easily  be  solved  by  iterations.  For  detona- 
tions, for  which  nrf   is  rather  small,  the  equation  can  in  good 
approximation  be  replaced  by 


2T 


(^+1)' 


^  '  %  *  2(^+1)  ^—  9^  "  ?Tl  9  ^  • 

Once    t,    is  found  one  obtains    t,    from  P^/p     «    I 
=  (1  +  ^07  -  2)?^)/(  Y+  1),   further   r^/^v    =    V^     •      Finally 
one  obtains  v^  and  v^  from  (4),     Thus  all  quantities   character- 
izing a  Chapman-Jouguet  reaction  are  determined. 
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Appendix  III 

Von  Neimanp's   argument  In   support  of 
the    Chapman-Joxiguet   condition 

A  significant  argument  In  order   to    justify  the   Chapman- 
Jouguet   condition  was  given  try  von  Neumann    [16]  •      It  Is 
essentially  cased  on  the  fact  (disregarded   so  far  In   the 
present  report)   that  the  reaction  actually  does  not  take   place 
In  a  sharp  discontinuity  front  but  rather  In  a   zone  of  finite 
uldth.     In  the   terms  employed  before,  the  result  of  this   argument 
can  be  formulated  as  follows:     Strong  deflagrations  and  weak 
detonations  do  not  exist  provided   that  the   reaction  process 
Is  strictly  one -dimensional  and   stationary  when  observed  fran 
the  moving  reaction  zone,    and  provided    that  a   certain  additional 
condition  Is   satisfied. 

To  see   this   let  us   consider  an  explosive   confined   to 
a  tube  extended  along  the  x-axls.      The  unexploded  gas   contained 
in  the   space  x  <  0  and  being  In   the   state    (0)  moves  with  the 
velocity  V     Into   ths  cross-section  x  -  0  at  which  a  deflagration 
reaction  Is   to  start.     While  in  the  preceding  parts  of  this 
report  the  transition  from  the  state    (0)   to  the   completely 
burnt  state    (1)  was   considered   a  discontinuous   process,   we 
now   turn   to  a  more  refined  discussion  visualizing  the   transition 
as  continuous,   gradually  leading  from  the  state    (0)   through 
Intermediate   states    In     which  the  gas   Is  Incompletely  burnt, 
to   the  final  state    (1).      The  reaction  zone  may  extend  from 
the   cross-section  x  -  0  to   the  cross-section  x  =  a  out  of  which 
the   completely   burnt  gas  moves  with  the  velocity  v^  Into   the 
space   X  >  a.      At  each  cross-section  of   the  reaction   zone, 
where,   as  assumed,   all  quantities  remain  stationary,    the  gas 
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is  a  mixture  of  a  certain  fraction,  n,   of  burnt  gas   and  a 
fraction,   1  -  n,   of  still  unburnt  gas.     We   assume   that  the 
fraction  n   increases  monotonically  on  passing  across   the 
reaction  zone  from  the  value  n  =  0  at  x  s=  0  to   the  value 
n  ^  1  at  X  f=  a.      The  value  of  n  then  completely  characterizes 
the   cross-section  in  the  reaction   zone. 

It  is  an  important  point  of   the  argument,   originally 
suggested  by  G.   I.   Taylor,    that  each  section  of  the  reaction 
zone  betveen   the  unburnt  state    (0)   and  the  intermediate   state 
(n)   is   subject  to  the   same   laws  of  conservation  ?ihich  so  far 
•■ere  analyzed  in  Section  1  for  a   complete  reaction  process. 
The  mixture   of  burnt  and  unburnt  gases   characterized  by  the 
fraction  n  possesses  a   certain  energy  function  E  =  S   {©) 
and  therefore   the   same    conclusions   concerning  possible   tran- 
sitions can  be  drawn  as  were  drawn  in  Section   1  for  a  completed 
reaction.      Therefore,   a  relation  between   the  possible  values 
of  pressure   p  and  specific  volume    r    in  the   state   (n)   can  be 
established,   and  accordingly  we  may  draw  a  family  of   transition 
curves   in  a   (T,p)-plane  as   in  Figure  21,  each  member  corre- 
sponding  to  a  value  of  n  between  0  and  1.      The  curve  for 
n  -  0  is   the   shock  transition  curve;     it  passes   through  the 
point  (t  >p   ).     The  curve  for  n  =  1  contains   the  values 
(T.,p.)   characterizing  the   state    (1).     We  now   introduce   the 
additional  condition,   announced  before*      the  family  of 
transition  curves  has  no  envelope   and  behaves  as  indicated 
in   the  diagram  Figure  21. 
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h/ea^  ae^^ppra^ i  o-n 


hapr>nsyi  -Jcup  u€  f    a/ef/sfr^tiOTo 


Figure  21 

Family  of  transition  curves  in  the  (T,p) -plane  for 
incomplete  deflagrations,  fractions,  n,  of  conpletion 
ranging  from  n  =  0  to  n  =  1. 


Since  for  all  these  transitions  the  initial  state 
(t  ,p  )  and  the  initial  velocity  v^  are  the  same,  formula 


(4) 


2_2 


is  applicable;      it  states   that  the   points    (t,p)  for  all 
intermediate   transitions   lie  on  a   fixed  straight  line,  A  • 
The   points  of  intersection  of  this   straight  line  *ith  the 
transition  curves    (n)  determine   the   actual  partial  transition 
from   (0)   to    (n).      On  the   line  A  we  follow   the   set  of  points 
of   intersection,    letting  n  range   from  0   to    Ij      thus  *e   obtain 
a  well-defined  point  of   intersection   on  the   curve  n  =  1.      (Se€ 
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Figure  21).  By  virtue  of  the  assumption  that  the  family  of 
transition  cu3?vea  behaves  as  is  indicated  in  the  diagram. 
Figure  21,  no  other  point  of  intersection  of  the  line  A  with 
tbB  curve  n  =  1  can  be  reached  through  intermediate  points 
on  the  line  A  ,  since  this  would  imply  passing  through  a 
section  of  this  line  not  corresponding  to  partial  transitions 
for  n  <  1.  Now  it  is  obvious  that  the  points  on  the  transition 
curves  n  =  1  that  can  be  reached  in  this  way  correspond  to 
weak  or  to  Chapman-Jouguet  deflagrations.   Thus  strong 
deflagrations  are  excluded* 

This  result  is  sufficient  to  exclude  also  the  possibilities 
of  weak  detonations  under  the  assumptions  made.  As  was  mentioned 
in  Section  5,  p.  29,  a  detonation  process  can  be  considered 
a  deflagration  process  coalescing  with  a  pre -compress ion  shock 
preceding  it.  It  is  immediately  seen  that,  in  particular,  a 
weak  detonation  process  corresponds  to  a  shock  followed  by  a 
strong  deflagration  process.  Exclusion  of  strong  deflagrations 
therefore  implies  exclusion  of  weak  detonations. 

Inasmuch  as  strong  detonations,  under  proper  boundary 
conditions,  are  also  excluded,  the  Chapman -Jouguet  detonations 
remain  as  the  only  possible  ones. 

In  conclusion  it  should,  however,  be  said  that  the 
assumptions  of  one -dimensional,  i.e.  of  strictly  confined 
flow,  is  never  completely  satisfied  in  reality  for  detonation 
processes  and  that  the  possibility  of  three-dimensional  weak 
detonations  is  therefore  not  completely  excluded  hy  the 
reasoning  presented. 
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